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Abstract 

The observation of Higgs decays into heavy neutrinos would be strong evidence for new 
physics associated to neutrino masses. In this work we propose a search for such decays 
within the Type I seesaw model in the few-GeV mass range via displaced vertices. 

Using 300 fb -1 of integrated luminosity, at 13 TeV, we explore the region of parameter space 
where such decays are measurable. We show that, after imposing pseudorapidity cuts, there 
still exists a region where the number of events is larger than 0(10). We also find that 
conventional triggers can greatly limit the sensitivity of our signal, so we display several 
relevant kinematical distributions which might aid in the optimization of a dedicated trigger 
selection. 


1 Introduction 

The Type I Seesaw mechanism Una is possibly the simplest extension of the Standard Model that 
can explain the smallness of neutrino masses. Even though most realizations of this mechanism 
invoke extra sterile neutrinos with Majorana masses too heavy to be probed, the possibility that 
these masses lie at the electroweak scale range is not excluded, and could actually be a more 
natural scenario. Such a case generically requires small neutrino Yukawa couplings, of similar size 
as those of the light charged leptons, however, this is not the only possibility. An approximate 
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U(1)l lepton number symmetry can be imposed to protect the smallness of neutrino masses, 
allowing for larger Yukawa couplings and heavy masses at the electroweak scale [5H7]. These 
models imply new free parameters that cannot all be fixed by the light neutrino mass matrix, it 
is therefore of upmost importance to search for complementary tests. 

When the new heavy neutrinos are lighter than the Higgs, the latter can present novel decay 
channels, in particular, a decay into a light and a heavy neutrino [8] . This would be followed by a 
subsequent decay of the heavy neutrino via a charged or neutral current interaction. In a number 
of recent references, the study of such Higgs decays at the LHC has been performed, focussing on 
the decay channels N —> £ + £~u |[9] and N —>■ tqq' [10], 

If the heavy neutrinos have masses of the order of a few GeV, the Higgs decay can lead to a 
noticeable displaced vertex, which is potentially a very powerful signal to look for m This mass 
range is particularly interesting, because it might lead to successful baryogenesis T2 ~HT1. 

Recently, the putative signal of a displaced vertex from heavy neutrinos produced in W decays 
has been studied p3ijT5]- In contrast to the latter work, in this paper we consider the signal 
of displaced vertices at the LHC, resulting from Higgs decays to heavy neutrinos. Such a mea¬ 
surement would allow us to directly probe the neutrino-Higgs coupling, giving a strong signal in 
favour of the Type-I seesaw model. This signal is obtained within the framework of a minimal 3+2 
neutrino model with an approximate U(1)l symmetry^, and after imposing all existing constraints 
from neutrino masses, direct searches, neutrinoless double beta-decay and lepton flavour violating 
processes involving fi ey e transitions. 

Although production of heavy neutrinos from Higgs decay is more limited statistically than 
that from W decays, the two are sensitive to different combinations of parameters in the seesaw 
scenario, and are therefore complementary. The putative observation of both signals would be 
an unprecedented probe of the low-scale seesaw scenario. Conversely, the non-observation would 
impose stringent constraints that may be essential to rule out an interesting range of seesaw scales. 

The paper is organised as follows. In Section [2] we introduce the model and impose all existing 
constraints in a convenient parametrization. In Section [3l we review the Higgs decays to heavy 
neutrinos and quantify the size of the corresponding branching ratios in the presently allowed 
parameter space. In Section [2 we consider Higgs production from gluon fusion at the LHC and 
study the displaced vertex signature. We illustrate the reach of an LHC search on the parameters 
space of the model, and discuss the impact of several kinematic cuts. In Section [5] we conclude. 

Some useful formulae are presented in the appendices: in Appendix [A] the contribution to 
neutrinoless double beta decay and in Appendix [B] the differential decay rate of the Higgs into a 
heavy and a light neutrino in the lab frame. 

2 Parametrization and Constraints 

A minimal 3+2 neutrino model is characterized by the addition of two heavy sterile neutrinos. 
This is translated into a 5 x 5 neutrino mass matrix, which for the normal hierachy can be written 
in diagonal form as follows: 

M v = U* diag(0, m 2 , m 3 , Mi, M 2 ) U ] , (2.1) 

1 Similar analyses in more complicated models can be found, for instance, in mm and Contribution 18 of [T8] . 
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where M. u is in the basis where Y e is diagonal, and rrii are mass ordered. The parametrization 
of [T9] decomposes U into four blocks: 


U§x5 ~ 


{Uat) 3x3 {Uah) 3x2 


{U s£ ) 2x3 {U s h) 2x2 

For the normal hierarchy, each block can be parametrized in the following way: 

1 0 \ ( 0 


U a l — TpMNS 


0 H 


U a h ~ i TpMNS 


V Hm] /2 WM~ lf 2 


Ust = i[ 0 HM~ 1 / 2 Rm ] 2 


U sh = H , 


where: 


H=(l + mJ 2 iT Mr 1 R - 1/2 


*h “ m £ 


- 1/2 


H = (I + M ~ 1/2 Rm £ R ] Mu 1/2 


-1/2 


( 2 . 2 ) 


(2.3) 


(2.4) 


In the previous equations, we have a unitary matrix TpmnS; which would correspond to the 
observed neutrino mixing matrix in the limit H —> I. The diagonal heavy (mostly sterile) neutrino 
mass matrix is denoted as M h = diag(M 1 , M 2 ). The other two light (mostly active) massive 
neutrinos have a diagonal mass matrix denoted by m £ = diag(m 2 , m 3 ) = diag( ^Am^, \J Am^). 
Finally, we have a complex orthogonal matrix R [20], which is parametrized as: 


( cos( 6>45 + ^745) sin(0 45 + i7 45 ) \ 
\ - sin( 0 4 5 + * 745 ) cos(0 45 + ^ 745 ) J 


Thus, the only free parameters left in the neutrino mass matrix are the angles $ 45 , 745 , the heavy 
neutrino masses M \, M 2 , and the two CP phases present in I/pmns- One can demonstrate that all 
of the mixing angles 6 t j can be restricted to the first quadrant. 

If we want to express our results for the inverse hierarchy, we need to re-write Eq. (12.11) taking 
into account the appropriate ordering. This leads us to a different mixing matrix, V, such that: 


M v = V* diag (m 2 , m 3 , 0 , M u M 2 ) T f . 


( 2 . 6 ) 


The reordering can be done through a permutation matrix acting on the active states. In blocks, 
we have: 


V a i — U pmns 
V s£ = i( HM ~ 1/2 


H 0 
0 1 


RmV 2 0 


Yah ~ i TpMNS 
V sh = H . 


H raj 1 V M h 1/2 


(2.7) 


Let us comment on the role of 0 45 and 745 . For large I 745 1 >2 — 3, the hyperbolic sine and 
cosine in Eq. (12.5ft give essentially the same result (modulo a sign), such that $45 behaves as an 
overall phase: 

Rh 45|»1 = ( ^ ^ ) cosh 745 e^« . (2.8) 
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Here, the ± refers to the sign of 745. As we can see, 6*45 can be factorised out of the mixing, and 
plays no significant role within the phenomenology of the model. Thus, the relevant parameters 
in this model are the two heavy masses Mi, M2, and the angle 745. Moreover, for fixed heavy 
masses, increasing 745 makes active-heavy mixing grow exponent ialhfE 

In order to understand better the active-heavy mixing in this limit, let us also assume H ~ /. 
In this case, for the normal hierarchy, we can write: 


Uu = (U ah ) n = ±Z? H ^- cosh 745 e^ 45 , 
Ugn = (U ah )i2 = i cosh 745 e^ 5 , 


where: 


Z £ — (£/pMNs)f 3 ± i\l -^-(UpMNs)e 2 ■ 


m 2 


( 2 . 9 ) 

( 2 . 10 ) 

( 2 . 11 ) 


This structure is similar to that found in the literature (see for instance [ 21 ]). In addition, up to 
corrections of O ( m^/Mj ), we can write the Dirac mass matrices, associated to the Yukawas, as: 


(m D ) n = I iZ/"r . W cosh745 .1 : ''' l|r ' , 
(mo)® = cosh745e ¥l(?45 . 


( 2 . 12 ) 

( 2 . 13 ) 


Thus, in this limit the size of the Dirac masses is exponentially enhanced with respect to the naive 
expectations of Seesaw models, mo ~ ^rn % M r As an example, for values of 745 ~ 7 , one would 
expect an enhancement of O (10 3 ). 

For the inverted hierarchy, the results are identical, but including: 


7 ih 


(I4>mns)^2 ± i 



(f^PMNs) 


11 


( 2 . 14 ) 


In the following, we shall not use any of the approximations above, and shall always take the 
exact form of U a h and V a h. Also, for definiteness, we set all neutrino oscillation parameters to 
their best-fit points as in [22] . and all CPV phases to zero. 


2.1 Constraints 

There are three relevant constraints on the parameter space explored in this work. These come 
from neutrinoless double beta decay (Ou/ 3 / 3 ), lepton flavour violation (LFV) and direct searches. 

Neutrinoless Double Beta Decay 

Currently, the strongest constraints on Ou/ 3/3 come from Germanium and Xenon experiments. 
On the Germanium front, the GERDA, HDM and IGEX experiments have combined their data, 
and determined a lower bound on the lifetime, T®J 2 > 3.0 x 10 25 yr. This corresponds to an 

2 Notice that the unitarity of the mixing matrix is kept under control by the H and H matrices. 
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Figure 1: Left: The contours show the maximum value of 745 allowed by the lack of observation 
of Right: Bounds from /1 — e conversion in nuclei, for the case of degenerate masses. The 

dark blue area is excluded by experiments with Au, while the light blue area can be probed in the 
future with A1 experiments. 


effective mass mpp < 0.2 — 0.4 eV [23] . On the other hand, for Xenon, KamLAND-Zen and 
EXO-200 have jointly imposed a lower bound of Ty 2 > 3.4 x 10 25 yr. This would correspond to 
mpp < 0.12 — 0.25 eV [24] , A list of future Oix/1/1 experiments can be found in [25j[26i|. 

The non-observation of 0 v/3/3 can put very strong limits on the active-heavy mixing. To 
calculate this observable, we use the formulae derived in Appendix [A] based on the work in [27]. 
The left panel of Figure [T] shows the maximum allowed value of 745 as a function of the two 
neutrino masses, for the normal hierachy, given the current bounds. 

Notice that the bounds vanish for degenerate neutrinos, as expected from the second term 
in Eq. flA.llIh It turns out that, for degenerate neutrinos, one can describe the neutrino mixing 
matrix with an inverse seesaw-like structure, making evident the existence of an approximate 
U{1)l symmetry. This symmetry constrains lepton-number violating processes from being too 
large. This was already discussed in [28] , Furthermore, in this limit, the light neutrino masses are 
protected from large loop corrections [29] , Thus, for the rest of this work, we shall consider the 
degenerate case, Mi = M 2 . 


Lepton Flavour Violation 

The most relevant processes constraining our parameter space are radiative LFV and /i — e con¬ 
version in nuclei. 

Radiative LFV processes include /i —>• ey and r —>■ £7 decays. For the former, the MEG 
experiment has placed an upper bound of BR(/i —> ey) < 5.7 x 1CT 13 [30], and the future upgrade 
expects to reach a value around 5 x 10 -14 . For tau decays, both the Belle and BaBar experiments 
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Figure 2: Constraints placed on the heavy neutrino parameter space, due to direct searches. We 
show constraints for \U e4 \ 2 , \U^ 4 \ 2 and \U r4 \ 2 on the left, center and right panels, respectively. 

have constrained their branching ratios. The strongest ones are given by BaBar, of BR(r —> ey) < 
3.3 x 10” 8 and BR(r —>- /iy) < 4.4 x 10~ 8 [31] • The future expected sensitivity for both channels 
at Belle If is of e>(l(r 9 ). 

Another important process is /i — e conversion in nuclei. The SINDRUM-I1 experiment has 
imposed limits on the conversion rate associated with Ti (4.3 x 10“ 12 ) [32], Au (7 x 10~ 13 ) [33] and 
Pb (4.6 x ICC 11 ) [M]. There exist several experiments which will attempt to probe lower values, 
such as Mu2e (Al, ~ 5.4 x ICC 17 ) [35], COMET (Al, ~3x ICC 17 ) [36], and PRISM/PRIME (Ti, 
e>(l(C 18 )) [37]. 

To calculate all these processes, we use the formulae of [38], and references within. We find 
that all observables give competitive constraints, but the most stringent, both now and in the 
future, comes from /i — e conversion. As none of these processes have yet been observed, an upper 
limit is imposed on y 45 , even for the degenerate case. This is shown on the right panel of Figured] 

Direct searches 

Finally, we also need to apply direct search bounds. Many experiments have tried to produce, 
and detect, these heavy neutrinos. Again, the lack of observation puts constraints on active-heavy 
mixing. Providing a faithful interpretation of each result on the 3+2 model is beyond the scope 
of this work. Thus, we shall take the bounds as reported in [35], shown in Figured] and apply 
them directly to our framework. 

The most important direct search constraints for this work are those of DELPffl [40j. One 
must note that, although all three bounds by DELPHI seem competitive, the importance of one 
or another shall depend on the predictions for the mixing within the model. For instance, for the 
normal hierarchy, one finds that generally \U e4 \ 2 is smaller than |f/^ 4 1 2 and [LC^I 2 by an order of 
magnitude, meaning that the latter two bounds shall be more stringent. 


3 Higgs Decays into Heavy Neutrinos 

As mentioned previously, observing Higgs decay into neutrinos would be a strong signal in favour 
of the seesaw model. The Higgs partial decay width into two neutrinos was initially calculated 
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in [8j, and can be written as: 


(3.1) 

where is the Higgs mass, and the scalar and pseudoscalar couplings are: 

S = -^ 2 - (("h'i + "H) ^[C'p]) 2 i p = (('"W ~ m »i) ^rn[C tJ }j 2 , (3.2) 

with Cij = J2k=i PkJPkj- Moreover, A (a, b, c) = a 2 + 5 2 + c 2 — 2a6 — 26c—2ac is a kinematic function, 
and oj = l/n\ for n identical final states. We find that the largest branching ratio happens for the 
decay into one light and one heavy neutrino: 

q 2 , on2 9 

T (h mNj) = f - fm h (1 - y 2 Y |^f , (3.3) 

oZTT TYl^ 

where ijj = Mj/rrih and, for the normal hierarchy: 

° ij = ' ( H 2 m] / 2 B)M ~ 1/2 ) ' ^ 3 ' 4 ^ 

For the inverted hierarchy, one shifts the (2 — 3) rows to the (1 — 2) rows. Here we see the 
very important fact that the PMNS matrix does not appear in the partial width. In particular, 
this means that our results shall not depend on the unknown Majorana nor Dirac CP phases. 
This is not the case for heavy neutrino searches involving the W boson, where the active-heavy 
mixing plays a central role due to the presence of charged leptons. This shows that, within this 
framework, the measurement of both h —> riiNj and W —> £Nj decays would be complementary, 
with the possibility to access the value of the CP phases. 

We plot the h riiNj branching ratio in Figure [3] for the normal hierarchy, along with 
constraints from LFV and direct searches. We find that branching ratios as large as 0.01 are 
generally ruled out by the former constraints. Moreover, if we want to work with branching ratios 
large enough to provide a signal at the LHC, the heavy masses cannot have values under a few 
GeV. 


T (h ViVj) = 


UJ 


8n mi 


-A 1/2 ( 


m 


hi 


5 1 


[m Vi + m v 


mr 


+ P 1 


(m h 


rri h 


mi 


4 Displaced Vertices from Higgs Decays 

The heavy neutrinos are not stable, and eventually decay through charged and neutral current 
interactions. The decay channels, width and lifetime can be found, for instance, in [39^31] . If the 
neutrino transverse decay length lies between 1 mm and 1 m, a displaced vertex signal could be 
recorded at ATLAS and CMS ji2bT7] . 

We now consider the possibility of observing such a signal, as a product of Higgs boson decays. 
The first step is to identify the region of interest, that is, one where the decay length is acceptable 
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Figure 3: Branching ratio for h —> nN. LFV 
and direct search constraints are shown in 
red and blue, respectively. Contours indicate 
branching ratios of ICC 2 , 10~ 4 , 10“ 6 and 10 -8 , 
from inner to outer curve. 
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Figure 4: Decay length, tjvC. Constraints are 
shown as in Figure [3j The region between 
dashed lines has 1mm < tjvc < 10 3 mm, dot¬ 
ted lines indicate 10~ 3 mm < tnc < 10 6 mm. 
The decay length decreases from left to right. 


and the Higgs branching ratio is not too small. The transverse decay length £n t is related to the 
heavy neutrino lifetime Tn through: 


Fv'/- 


\Pn t \ 

M, 


TnC , 


(4.1) 


where both the heavy neutrino mass Mj and transverse momentum p^ T are measured in GeV. 
Thus, for a given mass, the requirement of having a visible £n t puts constraints on \'Pn t \Tn■ 

In order to get an approximate idea of the region of interest, we plot in Figure Q] the decay 
length tnc. The shape of the curves can be understood by realizing that tnc oc M 3 |Cm | 2 and 
then taking the logarithm. One needs to be aware that these curves are given only to roughly 
illustrate the region where displaced vertices might be visible. The parameter which ultimately 
defines the region is the transverse decay length £n t , which depends on the transverse momentum 
with which the neutrino is generated. 

The Tnc lines are not parallel to the contour lines for the h —> riiNj branching ratio. This 
means that at some point the region of interest shall have a too small probability for h —> riiNj 
decay. Thus, we find that we require heavy neutrino masses between 2-20 GeV to be able to probe 
a displaced signature without significantly reducing the expected number of events. 

We now estimate the number of displaced vertices in this region. We take Higgs production 
through gluon fusion, gg —$■ h, followed by the decay h -P- riiNj. If we do not consider specific 
final states after heavy neutrino decay, nor any kinematical cuts, the event rate with a measurable 
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Figure 5: Region sensitive to events with a displaced vertex. The region in red would have more 
than 250 events with a displaced vertex, for an integrated luminosity of 300 fb” 1 , at 13 TeV. The 
region in orange would have more than 50 events. The blue region is ruled out by direct searches, 
and the dashed line indicates the reach of /i — e conversion experiments using A1 nuclei. The dot 
represents a benchmark point. 


displaced vertex at the LHC is: 


N = C I d \p hl , | dy h d \p Nr \ d(j) N 

d 2 a(gg —> h) d 2 T{h —y riiN, 


x 


d | p hT | dy h T h d \p Nx \ d(j> N 


—Q h (^v T ( mm ) - !) ( 1q3 - Cv T ( mm )) (4.2) 


Here, ph T and yh are the Higgs transverse momentum and rapidity. The azimuthal angle between 
Ph T and p Nr is denoted by In addition, T h is the Higgs width, while 7 h is the corresponding 
relativistic factor, 7/1 = Eh/rrih . The luminosity C is taken equal to 300 fb _1 for the second run of 
the LHC. The two Heaviside 0# functions make sure the decay length lies within the detection 
capability. The integration is performed with the Vegas subroutine of the CUBA library [48]. 

The integrand contains two differential distributions. The first one corresponds to gg —> h 
production, as a function of the Higgs transverse momentum and rapidity. This is obtained 
through the codes SusHi [49l[50lj and MoRe-SusHi [5T1I52] . The second is the Higgs differential 
decay width, in the lab frame, for decay into one light and one heavy neutrino. This is shown 
in Appendix [B] Note that we are using cylindrical coordinates. This is done in order to directly 
constrain the transverse decay length through the integration limits for \pn t \ ■ 

The region of parameter space leading to events with a visible displaced vertex is shown in 
Figure [5] for the normal hierarchy. We show the regions excluded by direct searches in blue, and 
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the reach of future LFV experiments by the dashed curve. The region in red would have more 
than 250 events with a displaced vertex, using 300 fb _1 of integrated luminosity at 13 TeV. The 
region in orange would have 50 events for the same luminosity. 

As mentioned previously, the left and right boundaries are determined by the experimental 
requirements on the decay length. For smaller masses, the decay length is too large, and the 
heavy neutrino escapes the detector. For larger masses, the decay length is too small, and the 
detector resolution is incapable of discriminating the displaced vertex from the interaction point. 
This constraint is imposed by introducing Eq. (14.11) . which depends on the heavy neutrino mass 
and mixing, within the Heaviside functions in Eq. (14.2[) . This result, also depends on the Higgs 
branching ratio being large enough, which determines the lower boundary of the region. 

Thus, we see that regions with heavy neutrino masses between 2-20 GeV favour displaced vertex 
events at the LHC, for values of |1 2 between O (10 _ ‘) and O (10~ 5 ). For future comparison, we 
establish a benchmark point for the normal hierarchy, with Mi = M 2 = 6 GeV and 745 = 8 . This 
point leads to 428 events with a detectable displaced vertex, and is displayed as a dot in Figure [5] 


4.1 Signatures from Heavy Neutrino Decays 

The previous section allowed us to calculate the number of events with a displaced vertex hap¬ 
pening due to Higgs decays. However, these events are not necessarily observable. The heavy 
neutrino eventually decays into other final state particles, which need to be detected. 

We obtain the differential decay rate for heavy neutrinos, and convolute it in Eq. (14.2[) . Since 
the heavy neutrino is lighter than the W boson, two-body decays are not allowed. Therefore, it 
will decay through a three-body process. Given its relatively large branching ratio, we focus on 
Nj —* pqq' decay, where the momenta of the final states is labelled by pi, p 2 and p 3 , respectively. 
To calculate the differential decay rate, we follow [53] . The procedure is carried out in two frames. 
First, part of the integration is done in the frame where the spatial component of pjy — p 3 vanishes. 
On this frame, the momentum components shall be denoted with a tilde (i.e. <pi). The rest of 
the integration is then performed in the frame where the heavy neutrino is at rest. Momentum 
components in this frame shall be denoted with a hat (i.e d 3 ). Finally, we relate these variables 
with the appropriate ones in the lab frame, as this is where the experimental cuts are placed. 
These shall be denoted with a “lab” superscript (i.e. p( ab ). 

The observed number of events is given by: 


N = C dYL 


d 2 a(gg —» h) 7 ^ d 2 Y(h —» n^Nj) 7 N dY(Nj —> /i qq') 

d\ph T \dyh Yh d\piy T \d(j)N Fjv dsd(cos 9^)d{cosdi)d^d^i 


cuts } 


(4.3) 


where the neutrino differential decay width can be found in [53] . and: 

dYl = d\ph T \ dyhd\pN T \d(j)Nds d(cos93)d(cos9i)d(p3d(j)i . (4.4) 

The integration is carried out over four new angular variables, as well as s = (pn — P3) 2 /M 2 . 

Furthermore, F cuts is a function describing all experimental cuts. For instance, if the experiment 
was to impose cuts on the transverse momenta and pseudorapidities of each final state, we would 
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Figure 6 : Similar to Figure |5j number of events in the h —> nN —>■ nqiqq' channel. The region in 
brown (green) would have more than 100 (10) events with a displaced vertex. The dash-dotted 
line indicates more than one event. On the left we show the region if no cuts are applied to the 
final states, on the right we only apply pseudorapidity cuts on all final states. 


have: 


Fcnts = Qh (^v T (mm) - 1) 0# (lO 3 -^jv T (mm)) (4.5) 

xn. i: 0 H (|p iT | lab - \p lT r n ) Gn(vr * - ?T) • (4-6) 

To get a better understanding of the sensitivity of this signal, we plot on the left panel of 
Figure [H] the region with visible displaced vertices, assuming no cuts on the final state. Our result 
essentially shows the same information as Figure [5l weighted by the N —> pqq' branching ratio for 
each point. This scales the number of events by a factor 1/4 - 1/5. The right panel of Figure [ 6 ] 
shows the same region, but including also conventional pseudorapidity cuts, that is, |r/ M | < 2.4 for 
the muon, and \rji\ < 2.5 for every other particle. We find that, although the overall shape of the 
region remains unchanged, the number of events is affected by the cut. 

As an example, we report the results for our benchmark point. In Figure | 6 l we have a total of 
110 events on the h —y riiNj —y qiqq' channel with no cuts (left panel), which is further reduced to 
78 events once all pseudorapidity constraints are applied (right panel). This is to be compared to 
the 428 events we expected from h —$■ riiNj decays (Figure EJ. 

4.2 Impact due to Kinematical Cuts 

For the purpose of giving a perspective of a future experimental search, we discuss the impact of 
several kinematic cuts on our analysis. 
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Figure 7: Ratio between the number of events with and without a cut on the muon transverse 
momentum, for the h riiNj —> riiHqq' channel. The blue (orange) points represent M t = 3 GeV 
(15 GeV). 


In the following, we plot the ratio of surviving events for each cut, imposing at the same time 
the displaced vertex and pseudorapidity constraints previously discussed. In order to understand 
the impact of the heavy neutrino mass, we show results for Mj = 3, 15 GeV, which are limiting 
values of our signal region. We find that the results are not strongly influenced by the heavy 
neutrino mass. For each cut, we also compare the exact number of events for our benchmark 
point. 

The first constraint we study is a cut on the transverse momentum of the muon, p™*;. This is 
shown in Figured We find that typical cuts between 20 and 30 GeV would reduce the number of 
events to a total between 40% and 20%. As an example, the benchmark point shows 32 (17) events 
after imposing a 20 (30) GeV cut. Since in this analysis we are not including detector effects, 
such as efficiency, it is clear that we need to relax the stringency of p™ 1 if we want to significantly 
improve the sensitivity with respect to that from DELPHI. We consider that, as these muons are 
not produced at the interaction point, a dedicated trigger with a smaller cut on the transverse 
momentum is more appropriate. 

Another cut of interest is that on missing transverse energy, p T = \p„ T \, shown on Figured 
We observe that one can impose $ T cuts up to 40 GeV without reducing the number of events 
below 80%. However, at this point the ratio drops, and one fold that for cuts above 70 GeV the 
ratio is again under 20%. For our benchmark point, imposing displaced vertex, pseudorapidity 
and $ T constraints, we find 75 (15) events for a $ T cut of 40 (70) GeV. 

Finally, in Figure d we show the impact of a cut on M eff = |p /iT | + |p gr | + \pq> \ + $t- As 
expected, we fold an endpoint for M eS = rrih . In this case it is possible to impose cuts as large 
as 95 GeV without reducing the number of events under 80%. Again, on our benchmark point, 
cutting on M e g = 95 GeV decreases the number of events to 63. 
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Figure 8 : As Figure [TJ but applying a cut on 
missing transverse energy. 


Figure 9: As Figure [TJ but applying a cut on 


5 Discussion 

In this work, we study the possible observation of Higgs decays involving heavy neutrinos, by 
means of a search for displaced vertices. 

This study is done in the context of the minimal 3+2 neutrino model, which is based on a 
Type-I Seesaw with two heavy sterile neutrinos. After imposing all constraints on the parameter 
space, we find that the model can be described in terms of two additional parameters, apart from 
the light neutrino masses and mixings. The two new parameters are a degenerate mass for the 
two heavy neutrinos, and the enhancement parameter in the Casas-Ibarra R matrix, 745 . 

We then calculate the partial width for Higgs decay into any two neutrinos. We find that 
the h —>■ n,Nj channel has the largest branching ratio, and concentrate on the description of a 
displaced vertex signal. This signal is particularly relevant for degenerate heavy neutrino masses 
of the order of a few GeV. 

It is important to stress that this prediction depends on the neutrino masses and on 745 , with 
no dependence on the neutrino mixing angles nor phases of the PMNS matrix. Such decays can 
therefore provide direct information on the new parameters of the model. 

For the LHC Run 2, there exist allowed regions of parameter space where the number of Higgs 
decays with a displaced vertex could be as large as 0 ( 100 ), before any other kinematical cut. 
It is important to note that the observation of the displaced vertex relies strongly on the decay 
channel of the heavy neutrinos and on the detection efficiency. As an example, we have included 
the branching ratio due to N —$■ /. iqq' decay, and imposed pseudorapidity cuts on the final states. 
Both considerations reduce the number of events down to O (20%) from the original number, still 
leaving a large enough amount to be observed. 

In order to perform a more realistic assessment of the signal strength, we have considered 
additional kinematical cuts. For instance, in the N —y (iqq' channel, we find that a 30 GeV cut on 
the muon transverse momentum, (typical of a level 1 trigger), the number of observable events is 
reduced to 4% of the initial number. This low efficiency is due to the low value of the momenta 
of final states, which in turn is a consequence of the low mass of the heavy neutrinos. Therefore, 
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as this does not include further potential losses from detector reconstruction inefficiencies, we 
conclude that one cannot rely on conventional cuts to properly observe this channel. In order to 
avoid this situation, we present two alternative kinematical cuts with much better efficiency. Such 
cuts are based on missing transverse energy, and effective mass, M e g. We believe this could 
be helpful in designing a dedicated trigger, and point out that such a trigger could be useful more 
generally to search for weakly interacting light particles. 

The measurement of Higgs decays to heavy neutrinos would constitute a powerful test of the 
mechanism of neutrino mass generation. This process can provide complementary information to 
the one that can be measured via the dominant production mechanism, W — y £Nj decays. 
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A Neutrinoless Double Beta Decay 

Following m, we find that the 0 v/3/3 amplitude A is proportional to: 

3 2 

^ V m « +E M - u k , (A.i) 

ti = 1 i= 1 

where Ai 0u/313 is the nuclear matrix element. Furthermore, for heavy neutrinos with mass larger 
than 100 MeV, one has: 

M 0 "W(m ei ) -> M°"W( 0) , (A.2) 

M 0u ^(Mi) -> 0 . (A.3) 

It is common in the literature to define: 

3 2 

mpp = m * U ^i = ~ Mi U *hi ’ A- 4 ) 

ti =i *=i 
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where the last equality is guaranteed by the seesaw mechanism, at tree level. Then, for non¬ 
degenerate masses, we can understand the heavy neutrino contribution by writing: 


2 

2—1 


f Y, M i u lhA - M 0 u W(M2)) M, U 2 hl 

-mpp M°^(M 2 ) + AM 0u ^M 1 U 2 hl (A.5) 


The hrst term on the right is proportional to the contribution from the light neutrinos, but 
strongly suppressed by the matrix element involving the heaviest neutrino. Thus, for heavy 
neutrino masses larger than 100 MeV, the second term would provide the dominant contribution. 
In this model, the matrix element U e h t can be written: 


U eh 


i 


Q e I qe r> t 

°sol n 'li ' °atm rL 2i 



For the normal hierarchy we have: 

(‘ S ' soOnH = [(^ PMNs ) l 2 - f^ll + ( h r pMNs ) l 3 hf 2 l ] ( Am ^) 1 / 4 
('S'atm)NH = [(^ PMNs ) l 2#12 + (^ PMNs ) 13 -^ 22 ] ( A '/ n ^ m ) 1/4 

while for the inverted hierarchy: 

(*-^sol)lH = [(^PMNs)llhfn + (t/pMNs) 12-^2l] (Am^) 1 / 4 
(^atm)IH = [(^PMNs)ll-f^l2 + (£A>MNs) 12 -^ 22 ] (Am atm) 

Then, the whole amplitude is proportional to: 

A oc 0) - Af°^(M 2 )) 

"A M 0 ^ [(S s e ol + 5£J + (S s e ol - S e at J e -<(®«-*r«)] 2 


(A. 6 ) 

(A.7) 

(A. 8 ) 

(A.9) 

(A. 10 ) 


(A.ll) 


Here, we see that the amplitude can be exponentially enhanced by 745 , negating the suppression 
from the matrix elements. For very large 745 , the only way to control this enhancement is by 
having degenerate heavy neutrino masses. 


B Higgs Decays into Heavy Neutrinos 

For completeness, we report the differential Higgs decay width, in cylindrical coordinates, on a 
boosted frame. We take a vanishing light neutrino mass, and take the heavy neutrino mass equal 
to Mj\ 


d 2 V(h ->• mNj) 
d \p Nr I d(fi N 


x 


1 \Pn t \ 

87f2 y/ m l + \Ph T \ 2 


(S + P)^ 


M[ 

m 2 h 


1 

1 

1 


Pn z {E 1 + E 2 ) 

~ Ph z Ei 

1 

PN Z =P Z 

Pn z (E 1 + E 2 ) — p hz E\ 

PN Z =Pz. 


(B.l) 
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Here, the energy of the outgoing neutrinos is: 


Ei = ^Mf + \p NT \ 2 +p 2 Nz (B.2) 

1 /2 

E 2 = [\pn t \ 2 +P 2 n z + \Ph T \ 2 +pI z -2\ph T \\pN T \ COS0JV -2ph z PN z ] ■ (B.3) 

The Higgs momentum in the direction of the beam axis is defined in terms of the transverse 
momentum and rapidity: 


Phz = \h n i + |iS. T | 2 Stall 2/k 


(B.4) 


Similarly, the variable p^ z is the heavy neutrino momentum on the direction of the beam axis. 
It is fixed by momentum conservation, and has the following allowed values: 


Pz = 


2 (™ 2 h + \Ph T | 2 ) 

± 


{m 2 h + M 2 + 21 p Nr 11 p hT | cos <; b N ) p hz 


« + \p hT \ 2 +P 2 hz )( (m 2 h - M 2 ) 2 - A(m 2 h \p NT \ 2 + M 2 \p hT \ 2 ) 


4 \pn t \ \Ph T I siir &N + 4(m^ + Mf )\pn t \ I Ph T I cos (, b N 


1/2 


(B.5) 


Demanding p% to have real values puts constraints on \pn t \ and (p^r- We find that, if: 

,2 n j2 \ 


(mi — M 2 ) 

\p hT \^ { h 


Alternatively, if: 


« ~ M 2 ) 

2 Af,- 


2 Mj 


< \Ph T \ 


o < \pn t \ < Pt 

—7r ^ (pN A tt 


Pt < |Pv T | < Pr 
sm to < TBTKTF 


(B.6) 


(B.7) 


The values are dehned as: 


Pt 


2(t7t| + |ph T | 2 sin 2 0 at) 


(m 2 h + M 2 )\p hl 


cos (p N 


± 


(m 2 h + \Ph T \ 2 ) x (("4 - M ]) 2 ~ 4M 2 \Ph T \ 2 sin 2 (p N ) 


1/2 


(B.8) 
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